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Abstract —This letter proposes a control theoretic approach to 
achieve proportional fairness amongst access categories (ACs) 
in an error-prone EDCA WLAN for provision of distinct QoS 
requirements and priority parameters. The approach adaptively 
adjusts the minimum contention window of each AC to derive 
the station attempt probability to its optimum which leads to 
a proportional fair allocation of station throughputs. Evaluation 
results demonstrate that the proposed control approach has high 
accuracy performance and fast convergence speed for general 
network scenarios. 


I. Introduction 

802.lie Enhanced Distributed Channel Access (EDCA) 
extends the basic Distributed Coordination Function (DCF) 
by classifying traffic into different Access Categories (ACs) 
to provide service differentiation. Traffic with higher QoS 
requirements, e.g. shorter delay deadline, is assigned a higher 
priority, and thus, on average, waits for less time before being 
sent to the channel. Although it provides QoS enhancement, 
EDCA is essentially unfair because it benefits high-priority 
traffic at the cost of low-priority traffic’s performance. If a 
network is loaded with a large proportion of high-priority 
flows, the channel will nearly be all occupied by high-priority 
flows, e.g. VoIP or video streaming flows, low-priority traffic 
however , such as email or web browsing data, will suffer se¬ 
vere starvation. Therefore, it is of significance to seek for a fair 
allocation of network resources (e.g. throughput, airtime and 
etc.) amongst different traffic types, and meanwhile, guarantee 
the specific QoS requirements and service differentiation. 

Recently m derives the proportional fair allocation amongst 
different ACs in an 802.1 le EDCA WLAN given distinct 
average delay deadlines, and proposes a control theoretic 
method to achieve it. In this letter we extend that work in a 
number of directions. Firstly, the analysis in C] assumes that 
transmissions are error-free while in this letter we relax this 
assumption. Packet losses are caused by either collisions or 
channel errors. Secondly, in (T| we fix CW max = CW m i n 
while in this letter due to the assumption of error-prone 
channels, packets need to be retransmitted and hence we set 
CWmax = 2 m • CWmin and m is the retry limit. Thirdly, 
the retransmission delay in this letter is not simply the col¬ 
lision delay as considered in (T), it should also include the 
delay caused by retransmission(s) of a TXOP burst which is 
determined by the packet error rate. 
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II. Proportional fairness in 802.1 1e EDCA WLANs 
A. Network model 

We consider a single-hop 802.1 le EDCA WLAN with one 
AP and n client stations. Traffic flows are classified into N 
ACs. We assume that each station only transmits traffic of 
a single AC, so no virtual collisions are considered in our 
analysis. The number of stations in the it h AC is n*. Stations 
in the it h AC are assumed to have the same packet error rate 
Pi, and transmit packets at the same PHY data rate r*. The 
packet size of all ACs is L bits. 


B. Station throughput 

We start with the analysis of station throughput under 
saturated network loads. The following parameters are defined 
for AC i: 

• CW^ in is the minimum contention window; 

• CW^ ax is the maximum contention window; 

• AIFSi is the duration of the Arbitration Inter-Frame 
Space, i.e. AIFSi = SIFS + ti x a, in which F is 
the number of included time slots, SIFS is the duration 
of the Short Inter-Frame Space and a is the duration of 
a physical time slot; 

• Mi is the number of packets contained in a TXOP burst. 

Due to the use of TXOP bursting, the RTS/CTS exchange 

mechanism is used to make fast recovery from collisions. If 
any packet within the TXOP burst fails due to channel errors, 
the burst is terminated and the station contends again for the 
channel to retransmit the failed packet. Therefore, a MAC 
time slot may be a PHY idle slot, a TXOP transmission or 
a RTS/CTS transmission. 

Let Ti denote the probability that a station in AC i attempts 
to transmit in a time slot, and r = {t*, i = 0, • • • , N — 1}. 
The station throughput of AC i is given by 

Sj(r) = 

PTE(Pld* xop ) 

N—l N—l 

P J a + £ ni Pj E(7f op ) + (1 - P 1 - £ mPf)T^i 

2=0 2=0 

( 1 ) 
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in which P 1 = n (i - nr is the probability of an idle 

2 = 0 


N-l 

time slot and Pj = r^( 1 — Ti) ni_1 Yl (1 — Tj) nj is the 

J=0,j/2 

probability that a station of AC i makes a TXOP transmission. 
T c °i = Trts + EIFS is the duration of a collision, where 
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EIFS is the duration of the Extended Inter-Frame Space, 
which is given by EIFS = T ACK +SIFS+DIFS. E{T* xop ) 
is the expected duration of a TXOP transmission from a station 
of AC i, 

HT? op ) = - Pi) k ~ X Pi ( k (-+ T °°) + T °i ) 

k=l \ r i J 

+ (1 - Pi) Mi (m* ( f - + T°°) + T^j 

where T°° = TpuYhdr + 2 SIFS + Tack is the protocol 
overheads associated with a single packet transmission within 
a TXOP burst, and T° = T rts + SIFS + T C ts + AIFSi 
is the protocol overheads associated with the transmission of 
the burst. E (Pldl xop ) is the expected payload transmitted in 
a TXOP burst of AC i, 

Mi-1 

E (Pldf op ) = Y, (! - Pi) k PikL + (1 - Pi) Mi MiL 
k=l 


C. Average delay 

Next, we will calculate the average packet transmission 
delay. As the delay is mainly dominated by the number of 
retransmission attempts caused by either collisions or erro¬ 
neous transmissions, we ignore the countdown and blocking 
delays considered in 0. The average delay of a TXOP burst 
consists of the following two expected delays: 

• Expected collision delay : Assuming that m is the retry 
limit and CW^ ax = 2 m ClkC n , the expected collision 
delay is 

D c°l =T col(j2j( P t C )\l- p f) + ( m +V(F) m+1 ) 

in which Pf is the conditional collision probability, 

N—l 

p? = i - (i - nC- 1 ft 0-~ T i) ni 

j=QJAi 

• Expected TXOP transmission delay : The expected TXOP 
transmission delay is calculated by multiplying the ex¬ 
pected number of transmissions by the expected duration 
of a TXOP burst. 

/ m+1 

D tx° P =E ( T te<>P) / ^ j(i - pcy (pfy- 1 (i - pf) 
v j= i 

+(m+i>(i 

in which PT = 1 — (1 — Pi) Mi is the probability of a 
TXOP failure due to channel noise. 

Combining the above two delays, the average delay of a single 
packet within a TXOP burst is therefore given by 

TJCOl jjtxop 

A = "mT- 1 - 1 - (2) 

E k ( l - Pi) k ~ x Pi + M i (i - Pi) Mi 

k= 1 


D. Relationship between Ti and CW 

The station throughput and average delay derived above 
are functions of station attempt probabilities. Next we will 
derive the relationship between t* and CW^ in based on the 
throughput model proposed in 0 but consider error-prone 
channels. The probability that a transmission fails due to either 
collisions or channel noise is given by 

Pf = l-(l-Pf)(l-P t E ) (3) 


The probability that the backoff counter is suspended due to 
a busy channel during the period of AIFSi is 


N—l 


pP = i — 


(1 - nr- 1 n (! - 


j =0 
jAi 
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in which t m ^ n is the minimum t value among all ACs. The 
relationship between r* and CW^ in is then given by 


1 . CWtnin (1 ~ P[) (1 ~ (2if ) m+1 ) \ 

2(1 - PP) 2 (1 - P?) (1 - 2 Pf) (1 - (Pf ) m + 1 ) J 

( 4 ) 


E. Proportional fair allocation 

In this section we aim at finding the optimal r to achieve 
proportional fair allocation of station throughputs given dis¬ 
tinct average delay constraints for different ACs. 

1) Utility function: The utility function is defined as the 
sum of the log of station throughputs 

N—l 

max U(t) :=yy ni log Sj(r) 
i =0 

s. t. Di(r) < di 0 < i < N — 1, 

0<Ti<l 0 < i < N — 1. 

in which di is the delay deadline for AC i. The station 
throughput is given by Eqn. ([!]), and the average delay is given 
by Eqn. ([2]). 

2) Solving the optimisation problem: We follow the algo¬ 
rithm described in (D to solve this optimisation problem. The 
steps are summarised as below: 

(i) Convert to a convex optimisation problem by making a 
log transformation r\i = log ; 

(ii) Introduce Lagrange multipliers using Karush-Kuhn- 
Tucker (KKT) conditions; 

(iii) Calculate the optimal Lagrange multipliers using a 
standard sub-gradient approach; 

(iv) Calculate the optimal r by plugging the optimal La¬ 
grange multipliers back to the Lagrangian. 


III. Centralised closed-loop control approach 

In this section, we design a centralized adaptive control 
approach to implement the desirable proportional fairness in 
real networks. Based upon the analysis in Section |II-E| the 
proportional fairness is achieved when the station attempt 
probability reaches its optimum value. The variable r is only 
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determined by CWmi n with CWm ax , AIFS and TXOP 
taking the recommended values. Our approach uses a multi- 
variable closed-loop control system to tune CWmin = 
{CWmin, ^ ® 0, • * • , N — 1} so as to drive the station attempt 
probability to its optimum. As the station attempt probability is 
hard to measure in real networks, we measure the transmission 
failure probability p F (r) = {pf(r), i = 0, ••• ,7V — 1} 
instead of r. The observed failure probability can be estimated 
as pf = „ p %i > i n which N} represents the number of 
packets of AC i that are received by the AP within a beacon 
interval with the retry bit equal to 1, and N® represents the 
number of received packets with the retry bit set to 0. The 
estimation is based on the assumption that all packets will be 
eventually received within the retry limit. 

A. Adaptive control algorithm 

The adaptive control algorithm is performed centrally at the 
AP every beacon interval, which is typically 100ms in 802.11 
WLANs, and consists of two steps: 

1) At the end of each beacon interval, the observed failure 
probability p F resulting from the current CWmin is 
measured by the AP. 

2) The controller computes a new set of CWmin based 
on the measured p F . An update of CWmin is then 
distributed to the WLAN by the next beacon frame. 

B. Linearisation of the non-linear plant 

As the period of the adaptive control algorithm is long 
enough to assume that the measurement corresponds to sta¬ 
tionary conditions, i.e. the system has no memory, p F depends 
only on the current CWmin. Eqn. ^ and Eqn. 0 give a non¬ 
linear relationship between p F and CWmin • To simplify the 
controller design, we will use a linear approximation to this 
non-linear relationship around the stable point of operation, 

P F - P F * = ( cw min - CW^ in ) • H (5) 

in which 



d Po 

dp? 

9 Pn- i 

9CW^ in 

£>Po 

dCW fin 

dp f 

9CW min 

9 Pn- i 



9CW rnin 

9CW min. 


9 Po 

dp f 

9 Pn -i 

dCW f~n 

9CW*" 1 

9CWZ- 1 


C. State feedback control 

With the linearisation, the WLAN can be represented as a 
discrete MIMO state-space model. The failure probability at 
instant fc+1 is determined by the minimum contention window 
input to the WLAN at instant k, the state and measurement 
equations are therefore given by 

f x(k + 1) = Bu(k ) 

\ y{k) = Cx(k) 

in which the system state is x(k) = [_ p F (k )] T ; the system 
input is u(k) = [CWmin(k)] T ; the system model matrices 
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are B = H r and C = Inxn ; the system output is y{k) = 
Cx(k ) = [p F (k )] T . 

The control task can be accomplished by using the LQI 
control method |3) to design our controller. Fig. [T] shows the 
control block diagram for the system, in which the controller 
input r(k) = p F is the optimal failure probability, e(k) = 
r(k) — y(k) is the error between the output feedback and the 
input reference value, and s(k) = s(k — 1) -\-T s e(k — 1) is the 
output of a discrete integrator and T s is the sampling period. 
The K G M 7Vx2Ar is the control gain matrix. 

The LQI controller computes an optimal state-feedback 
control law 

u(k) = —Kz(k) 

by minimising the quadratic cost function 

oo 

j( U (k)) = J2 (z T (k)Qz(k) + u T (k)Ru(kj) 
k =o 

for any initial state as(0), in which z(k) = [x(k); s(k)]. The 
optimal state feedback gain matrix K is computed by solving 
the associated discrete algebraic Riccati equation El- 

The matrices Q and R are the weighting matrices respec¬ 
tively indicating the state and control cost penalties. Q and 
R are required to be real symmetric and positive definite. 
The selection of weighting matrices Q and R affect the 
convergence speed of the controller. 


IV. Performance evaluation 

A. Impact of packet error rates 

We will first evaluate the impact of packet error rates on 
the station throughput and delay performances. Consider a 
WLAN with four stations and each station transmits a traffic 
flow of a distinct AC. The protocol parameter values used 
in the simulations are listed in Table [I] Fig. [2] plots the 
throughput, average delay and station attempt probability of 
four ACs as the packet error rate pvi varies while keeping 
Pbe = Pvo = Pbk = 10 -6 . It can be seen that as the 
PER increases, the station in AC_VI attempts harder with an 
increasing attempt probability but the throughput and delay 
performances irreversibly degrade. The increase of pvi also 
results in the decrease of attempt probabilities of the other 
three ACs, but the coupling of station transmissions lead to 
invariant delay and throughput performances, which implies 
that the delay and throughput of an AC are only dependent on 
its own PER. 
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TABLE I 

802.11 PROTOCOL PARAMETERS USED IN THE SIMULATIONS 



Racket error rate of ACVI 


Fig. 2. Delay and throughput performances, L = 8000 bits, tlbe = n vi — 
nvo = n BK = 1, r B E = ryi = r vo = r B K = 54Mbps, p B E = 
Pvo — Pbk — 10 — cLbe — 400ms, dyj = 250ms, dyo — 200ms and 
dsK = 400ms. The delay deadlines are not reached in this example. 



Fig. 3. Air-time allocation, L = 8000 bits, tibe = nvt = nvo = 
n B K = 1, r B E = r VI = r vo = r bk = 54Mbps, pbe = Pvi = 
Pvo — Pbk — 0.001, dBE — 400ms, dyj = 210ms, dBK — 400ms. 
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Fig. 4. Injection and/or removal of stations in the WLAN 


B. Air-time 

The flow total air-time is defined as the fraction of time 
used for transmitting a flow. For a flow of AC i, the flow 
total air-time consists of the successful flow air-time and the 
collision air-time, given by 


It can be seen that CW m i n converges to the desirable value 
very quickly when the network size changes as long as proper 
Q and R are chosen. Moreover, the steady-state errors can be 
neglected, which means the control system has high accuracy 
performance. 


pair _ 


PTe(T- xop ) + Ti Pf ol T co ‘ 


pidleg- . 


N—l 

■ E mP^m 

i =0 


itxop 
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Again we consider a WLAN with four stations individually 
carrying a flow of a distinct AC. Fig. |3(a)| plots the flow total 
air-times of four ACs as the delay deadline of AC_VO varies 
while keeping that of the other three ACs fixed. It can be seen 
that when dyo < 210ms, the delay constraint of AC_VO 
keeps tight, and four flows have different air-times. As dyo 
increases, the air-time allocation is getting more and more 
even. When dyo increases up to 210ms, the four flow total 
air-times are equalised and all the delay constraints become 
loose now. Fig. |3(b)| plots the sum of flow total air-times. We 
find that if any of the imposed delay constraints is tight, the 
air-times sum is larger than 1, otherwise they sum to unity. As 
the delay deadline is increased, the sum of collision air-times is 
decreasing while the sum of successful air-times is increasing. 
The fact that the air-times sum is not less than 1 is reasonable 
as the flow air-time usage overlaps due to collisions. This does 
not imply that the channel idle probability P ldle = 0. 

C. Adaptivity 


..ACVO 
- ACBK 


-AC.BE 
- ACVI 
■ i ACVO 
.. AC.BK 


(a) Contention window 


(b) Throughput 


Fig. 5. Adaptivity to changes of network size, L = 8000 bits, pbe = 
pbk = 36Mbps, ryi = ryo = 54Mbps, pbe = Pbk = 0.001, pyi 
pyo = 0.01, dBE — dyj = dyo = dsK = 400ms. 


V. Conclusions 

This letter proposes a closed-loop control approach to 
achieve proportional fair allocation in an error-prone EDCA 
WLAN. The derived proportional fairness assigns unequal air¬ 
times to flows if the imposed delay constraint is tight. We have 
demonstrated using evaluation results that the proposed control 
approach has high accuracy and fast convergence speed, and 
is adaptive to general network scenarios. 


We will next evaluate the adaptivity of the proposed control 
method to changes in the network size. The example is 
depicted in Fig. [4] Four stations, one in each AC, join the 
network at t = 0s. The station in AC_VO leaves at t = 60s. 
One more station in AC_BE joins the network at t = 110s, 
and the station in AC_VI leaves at t = 210s. Fig. [5] plots the 
variation of CW m i n and the station throughput over time. The 
selection of matrices Q and R uses the trial and error method. 
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